On Asymptotic Stability of Solitary Waves 
in a Nonlinear Schrodinger Equation 



V. S. BuslaevS 

Department of Mathematical and Computational Physics 
Faculty of Physics 
St- Petersburg University, Russia 
e-mail: buslaev@mph.phys.spbu.ru 

A. I. Komechd 

Fakultdt fiir Mathematik, Universitdt Wien 
Nordbergstrasse 15, Wien A- 1090, Austria 
e-mail: alexander.komech@univie.ac.at 

E. A. Kopyloval 

Keldysh Institute for Applied Mathematics RAS 
Miusskaya sq. 4-, Moscow 125047, Russia 
e-mail: ek@vpti.vladimir.ru 

D. Stuart 

Centre for Mathematical Sciences, 
Wilberforce Road, Cambridge, CBS OWA 
e-mail: D.M.A.Stuart@damtp.cam.ac.uk 



Abstract 

The long-time asymptotics is analyzed for finite energy solutions of the ID Schrodinger 
equation coupled to a nonlinear oscillator. The coupled system is invariant with respect 
to the phase rotation group U{1). For initial states close to a solitary wave, the solution 
converges to a sum of another solitary wave and dispersive wave which is a solution to the 
free Schrodinger equation. The proofs use the strategy of Buslaev-Perelman [Sid]: the 
linerization of the dynamics on the solitary manifold, the symplectic orthogonal projection 
and method of major ants. 
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1 Introduction 



Our main goal is the study of the distinguished dynamical role of the "quantum stationary 
states" for a model ?7(l)-invariant nonlinear Schrodinger equation 

t) = -tp"{x, t) - 5{x)F{ip{0, t)), xe IR. (1.1) 

Here ip{x,t) is a continuous complex- valued wave function and F is a continuous function, the 
dots stand for the derivatives in t and the primes in x. All derivatives and the equation are 
understood in the distribution sense. Physically, equation (11.11) describes the system of the free 
Schrodinger equation coupled to an oscillator attached at the point x = 0: F is a nonlinear 
"oscillator force" . 

We assume that F{ip) = —VUlip) where U{ip) = u{\ip\). Then (II. ip defines a f/(l)-invariant 
Hamilton system and admits finite energy solutions of type ^/'^(x)e*'^* called solitary waves or 
nonlinear eigenfunctions. The solitary waves constitute a two-dimensional solitary manifold in 
the Hilbert phase space of finite energy states of the system. We prove the asymptotics of type 

7/>(-, t) ~ V'-±e^"±* + iy(t)$±, t ±oo, (1.2) 

where W{t) is the dynamical group of the free Schrodinger equation, ^± G Cb(lR) fl L^(IR) 
are the corresponding asymptotic scattering states, and the remainder converges to zero as 
C(|t|~^/^) in global norm of Cb(lR) fl L^(IR). Here Cb{]R) is the space of bounded continuous 
functions IR — C The asymptotics hold for the solutions with initial states close to the stable 
part of the solitary manifold, extending the results of [Sj HJ [32l |3ll [35] to the equation (11.11) . 

Let us note that we impose conditions which are more general than the standard ones in 
the following respects: 

i) We do not hypothesize any spectral properties of the hnearized equation, and do not require 
any smallness condition on the initial state (only closeness to the solitary manifold). 

ii) The stable part of the solitary manifold is characterized by a condition on the nonlinearity 

(13.41) . The relation of this to the standard criterion for orbital stability di^J \-ipi^{x)\'^dx > (see 
|11] and references therein) will be discussed below. 

This progress is possible on account of the simplicity of our model which allows an exact analysis 
of all spectral properties of the linearization. 

Our investigation is inspired by a fundamental problem of quantum mechanics. The solitary 
waves were introduced by Schrodinger for the quantum electron coupled to the Maxwell field, 
[38]. He identified the solitary waves with the quantum stationary states. The asymptotics 
of type (II. 2p for the coupled f/(l)-invariant Maxwell- Schrodinger equations would describe the 
transitions between the quantum stationary states, while the dispersive wave W{t)^± would 
correspond to the electromagnetic radiation. Let us note however, that the asymptotics of type 
(ll.2p are not proved yet for the coupled equations. 

For the first time, the asymptotics of type (II. 2p were established by Soffer and Weinstein 
[391 HO] (see also [36]) for nonlinear f/(l)-invariant Schrodinger equation with small initial states 
if the nonlinear coupling constant is sufficiently small. The next result was obtained by Buslaev 
and Perelman [3, ll| who proved that the solitary manifold attracts finite energy solutions of 
a ID nonlinear f/(l)-invariant translation invariant Schrodinger equation with initial states 
sufficiently close to the stable part of the solitary manifold. 

The novel techniques of Buslaev and Perelman are based on the symplectic geometry in 
Hilbert space and the spectral theory of nonself adjoint operators. These techniques were moti- 
vated by the investigation of soliton asymtotics for integrable equations (a survey can be found 
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in [7] and [10]), and by the methods introduced by Soffer and Weinstein [391 HOI El- Similar 
techniques were developed by Miller, Pego and Weinstein for the ID modified KdV and RLW 
equations, [32l [311 ES]- The methods and results were extended in [5l [6] to the Schrodinger 
equations with more complicated spectral properties, and in [T7] to a translation invariant 
wave-particle system. Further references can be found in [0] and [T7J. 

Let us comment on the general strategy of our proofs. We develop the approach [H [T7] for 
our problem. Firstly, we apply the symplectic projection onto the solitary manifold to sepa- 
rate the motion along the solitary manifold and in transversal direction. Secondly, we derive 
the modulation equations for the parameters of the symplectic projection, and linearize the 
transversal dynamics at the projection of the trajectory. The linearized equation is nonau- 
tonomous, and this is one of the fundamental difficulties in the proof. This difficulty is handled 
by the introduction of an autonomous equation (by freezing the time) with an application of 
the modulation equations to estimate the resulting additional error terms. A principal role in 
the rest of the proof is played by the uniform decay of the frozen linearized dynamics projected 
onto the continuous spectrum, and the method of majorants. 

Let us note the following two main novelties in our approach to the uniform decay. First, we 
calculate exactly all needed spectral properties of corresponding generator. Second, we do not 
use a spectral representation of the generator. Instead, we develop the Jensen-Kato approach 
applying directly the Zygmund type Lemma 6.1 (cf. [18\ Lemma 10.2]) to the Laplace integral of 
the resolvent. We expect that the development would be promising for more general problems. 

The paper is organized as follows. In Section [2] some notation and definitions are given. 
In Section [3] we describe all nonzero solitary waves and formulate the main theorem. The 
linearization at a solitary wave is carried out in Section HI In Sections |5] and [6l we construct the 
spectral representation for the linearized equation. In Section [7] we establish the time decay for 
the linearized equation in the continuous spectrum. In Section M the modulation equations for 
the parameters of the soliton are displayed. The decay of the transverse component is proved 
in Sections fTOfTTl In Section [12] we obtain the soliton asymptotics (11.21) . In Appendix we study 
the resolvent of linearized equation. 

In conclusion, we expect that the asymptotics (11.21) holds for any finite energy solution of the 
equation (11.11) . however this is still open problem. We hope to prove it combining our methods 
with the techniques of the papers [25| |2S], where global attraction to the solitary manifold is 
proved for the ID Klein-Gordon equation with the same delta-nonlinearity. We also intend 
to treat in a later publication the case when the linearization has nontrivial stable oscillatory 
modes (which occurs if (13.41) holds but (13. 5p does not). 

2 Notation and definitions 

We identify a complex number tp = tpi +itp2 with the real two-dimensional vector {ipi,ip2) £ IR^ 
and assume that the vector version F of the oscillator force F admits a real-valued potential, 

F{^) = -VU{^), U e C2(IR2). (2.1) 

Then (II. ip is formally a Hamiltonian system with Hamiltonian 

W) = ^/ \ijTdx + U{m)- (2.2) 
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which is conserved for sufficiently regular finite energy solutions. We assume that the potential 
[/('?/') satisfies the inequality 

U{z) >A- B\z\^ with some A e IR, 5 > 0. (2.3) 

Our key assumption concerns the ?7(l)-invariance of the oscillator (cf [2]), where U{1) stands 
for the group e*^, 6 G [0, 27r] acting by phase rotation ip ^ e^^ip. Namely, we assume that 
U{'^) = MdV'H with u G C2(1R). Therefore, by (EI]), 

F(^) = adV^n^;, V^gC, aGC^(lR), (2.4) 

where a{\il}\^) is real. Then F{e^^il}) = e^^F{%l)), 9 G [0, 27r] and F(0) = for continuous F. 
Obviously, the symmetry holds true if U{%1)) = ■u(|'?/'p). The symmetry implies that e^^ip{x,t) 
is a solution to (11.11) if ^J^x^t) is. The equation is [/(l)-invariant in the sense of [H], and the 
Nother theorem implies the charge conservation: 

Q(^) = j \^\^dx = const. (2.5) 

The main subject of this paper is an analysis of the special role played by "quantum stationary 
states" , or solitary waves in the sense of [H] , which are finite energy solutions of the form 

^(x,t) = ip^{x)e''^\ CO eM. (2.6) 

The frequency u and the amplitude ipuj{x) solve the following nonlinear eigenvalue problem: 

- uM^) = -^'^{x) - 5{x)F{^U0)), xeJR. (2.7) 

which follows directly from (II. ip and (12.41) since G IR. 

Definition 2.1 iS denotes the set of all nonzero solutions G H^{JR) to l[2. 7| j with all 

possible G IR. 

Here if^(lR) = denotes the Sobolev space of complex valued measurable functions with 
j {\ip'\'^ + \ip\^)dx < oo. We give below in section [3] a complete analysis of the set S of all 

nonzero solitary waves by an explicit calculation: it consists of functions C{uj)e^^^^^^^^ 

with C > 0, u; = uj{C) > and any 9 G [0, 27r], and C restricted to lie in a set which, in the case 
of polynomial F, is a finite union of one-dimensional intervals. Notice that C = corresponds 
to the zero function iIj{x) = which is always a solitary wave as F(0) = 0, and for < only 
the zero solitary wave exists. 

Our main results describe the large time behavior of the global solutions whose existence is 
guaranteed by the following theorem, which is proved in [24J. 

Theorem 2.2 i) Let conditions h2. h2.3\) and ^2.4^ hold. Then for any ipo{x) G there 
exist a unique solution ip{x, t) G C;,(1R, H^) to the equation (QHP with initial condition ip{x, 0) = 

ii) There exists A{tpo) > such that the following a priori bound holds: 

sup||V^(t)||^i<A(V^o)<oo. (2.8) 

The functional spaces we are going to consider are the weighted Banach spaces L^, p G 
[1, oo), /5 G IR of complex valued measurable functions with the norm 

II^IU. = ||(1 + |x|)M^)IIl.. (2.9) 
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3 Solitary waves and statement of the main theorem 

Lemma 3.1 The set of all nonzero solitary waves is given by 

S = {^^e'^ = Ce^^-v^l^l : uj>0, C > 0, ^= a{C^)/2 > 0, de [0, 27r]}. 

Proof Let us calculate all solitary waves (12.61) . The equation (12. 7p implies ip"{x) = uip^x), 
X ^ 0, hence the formula ■ijj{x) = C±e^^ gives two linearly independent solutions in each of 
the two regions ±x > depending on which branch of ^/lj is chosen. Since iplx) G it is 
necessary that uj > and the branch is chosen with ±y/uj > for ±x < 0. Furthermore, since 
ip'{x) G L^, the function ip{x) is continuous, hence C_ = C+ = C and the solutions are of the 
form 

^/,(a;) = C'e-'^l^l, K = v^>0, iu>0. (3.1) 

Finally we get an algebraic equation for the constant C equating the coefficients of 5{x) in both 
sides of flO) : 

= ^'(0+)-^'(0-)+F(^(0)). (3.2) 
This implies = —2k,C + F{C), or equivalently. 



Corollary 3.2 The set S is a smooth manifold with co-ordinates 9 G lRmod27r and C > 
such that a^C^) > 0. 

Remark 3.3 We will analyse only the solitary waves with a'{C) ^ 0. On the manifold S we 
have uo = n'^ with k = a(C^)/2 according to liS. 3\) . Hence, the parameters 6,uj locally also are 
smooth coordinates on S at the points with a' = a'{C) ^ since ui' = 2nn' = aa'C ^ then. 



The soliton solution is a trajectory ipuj(t){x)e'^^^^'' = Ce V'^WI^Ie*^^*), where the parameters satisfy 
the equation 9 = uj, uj = 0. The solitary waves e*^'?/'^(x) map out in time an orbit of the U{1) 
symmetry group. This group acts on the phase space H^(R) preserving Ti. and the standard 
symplectic form (16.31) : in other words the solitary waves (12.61) are relative equilibria of the 
corresponding Hamiltonian system. 

Let us denote N{C) = J \'ijjujix)\'^dx with uj = k^, and k = a(C^)/2 according to (13.31) . It is 

easy to compute that N{C) = C^/k. We now differentiate: 

N'iO = ^ _ 

Differentiating the identity (13.31) . we obtain k' = a'C. Thus, again by (13. 3p . 

N\C) = ^(1 - ^) ^ 
K a 

if C > 0, a > and a' ^ a/C^. Therefore noticing that N'{C) = u'{C)d^ j \il)u^\'^dx with 
uj'{C) = 2k,k,' = aa'C, we obtain the following result 
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Lemma 3.4 For C > 0, a > we have 

J \ilJu;{x)\'^dx < z/ a' G (-00,0) U (a/C^+oo), 

and 

I \ij^{x)\'^dx > ifO<a'<a/C^. 

Remark 3.5 (i) Orbital stability of solitary waves is a much studied subject (see [H] for 
very general theorems in this area, and [13] for an approach more similar to that taken in 
this paper). The standard condition for orbital stability (|Il]) for the present problem would 

read (9^ J \ip^j{x)\'^dx > 0; this is expected to be a necessary and sufficient condition for orbital 

stability when the Hessian of the augmented Hamiltonian ([13]) has a single negative eigenvalue. 
In the present problem it can be easily calculated that this Hessian is non-negative when a' < 
and thus the standard condition is not necessarily relevant if a' < 0. Indeed Theorem 13.71 
asserts stability in the case a' < 0. Restricting to a' > 0, in which case the Hessian does have 
a single negative eigenvalue, the calculation above shows that orbital stability is expected to 
hold when a' < a/C"^. In this paper we will work under the spectral condition (13.51) which, for 
a' > 0, is slightly stricter: it is imposed to ensure that the linearization has no discrete spectrum 
except zero (which is always present on account of the circular symmetry of the problem). If 
a/\/2C'^ < a' < a/C'^ there are two purely imaginary eigenvalues of the linearized operator. 
It is intended to treat this case in a later publication thus extending our proof of asymptotic 
stability to the entire range 

- oo < a' < a/C^. (3.4) 

For a' > a/C"^ the linearized operator has a positive eigenvalue and the solitary wave is linearly 
unstable. 

(ii) It is explained at the end of section H] that (13.41) can be interpreted as saying the 
restriction of the standard symplectic form (16. 3p to the tangent space to S is non-degenerate 
(i.e. S satisfies the condition to be a symplectic submanifold) . 

Definition 3.6 We say the solitary wave ip^{x)e^^ = (7e^^l^l+'^, C > satisfies the spectral 
condition if u > and (cf. Remark l37^) 

a'{C^) e (-00, 0) U (0, a{C^)/{V2C^)). (3.5) 

Let us denote by W{t) the dynamical group of the free Schrodinger equation: W{t)f is de- 
fined by the Fourier representation for all tempered distributions /. Our main theorem is the 
following: 

Theorem 3.7 Let conditions K2l]} . and (2^ hold, P>2 andip{x,t) e C(IR,H^) be the 

solution to the equation U.l\) with initial value ipoi^) = ^'(^^jO) G fl Lj^ which is close to a 
solitary wave ipujoe'^° = CqC'v^I^'I+^^o with Cq > and iUo> 0: 

d:=Uo-^.oe'''\\mnLl'^l- (3.6) 

Assume further that the spectral condition 5\) holds for the solitary wave with C = Cq. Then 
for d sufficiently small the solution admits the following asymptotics: 

^(■, t) = V^.^e'^i* + + r±(t), t ±oo, (3.7) 

where $± G Cb(lR) nL^(lR) are the corresponding asymptotic scattering states, and 

\\r±{t)\\c,mnL^m = 0{\t\-^'^), t ^ ±oo. (3.8) 
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Remark 3.8 It is possible to derive further information about the structure of ^± and r±{t) 
as discussed towards the end of section 10. 

4 Linearization on the solitary wave 

As the first step in the proof of main theorem, let us linearize the nonhnear Schrodinger equation 
( 11.11) at a sohtary wave e*^'^*+^Va;(a;), with V'a;(a^) = Ce"'^'^' where k = y/cu > and C > 0. 
Substituting 

^(a;,t) = e*(-*+^)(^^(x)+x(x,t)) (4.1) 

to (11.11) . we obtain, 

- Luxix, t) + zxix, t) = -x"{x, t) - Six) [F{C + x(0, t)) - F{C)] (4.2) 
Use the representation (12. 4p to write 

FiC + x)-FiC) = a{\C + x\')iC + x)-ai\C\')C 

= a{{C + x)(C + mC + x)-a{\Cf)C 

= a{\C\')x + a\\C\')C{Cx + Cx) + 0{\x\') 

= a{C')x + a'{C')C\x + x) + 0{\xr) (4.3) 
since C > 0. Hence, the first order part of (14.21) is given by 
ix{x,t) = -x"{x,t) +ujx{x,t) 

~5ixMC')xiO,t) + a'iC')C'2RexiO,t)]. (4.4) 

Now it is evident that the first order part is not hnear over the complex field. On the other 
hand, it is linear over the real field. Hence, it would be useful to rewrite (14. 4p in the real form. 
Namely, identify X = Xi + ^X2 £ Cwith the real vector (xi, X2) £ ^ and denote it again by x- 
Then fl4.4p becomes 

3X{.x,t) = -x"ix,t) +ujxix,t) 



-6{x)[a{C^) + 2a'{C^)C^Pi]x{0,t), (4.5) 
X2 I I 

(4.6) 



where is the projector in IR^ acting as ^ PfM and ?' is the 2x2 matrix 



-1 

1 



Respectively, we also rewrite (II. ip in the real form 

jij{x, t) = -<(x, t) - <5(a;)F(^(0, t)), (4.7) 
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as an equation for ^(x,t) G IR^ with F{ip) G IR^ which is the real vector version of F(^) G C 
Then the hnearization (14 .Sp reads as 



jxix, t) = -x"{x, t) + iuxix, t) - 5(x)F'((C, 0))x(0, t). (4.8) 

Here F' is the differential of the map F : 1R2 ^ IRl 

F'((C,0)) = a + 6Pi, a:=a(C2), b := 2a'{C^)C^. (4.9) 

In order to apply the Laplace transform the next step is to complexify the system (14.81) 
i.e. to consider it as a system of equations for the complex functions Xi{x , t) , X2{x , t) , so 
x{x,t) G for any fixed {x,t). This gives a system which is linear over the complex field 
allowing application of the Laplace transform. To write this system more concisely let us 
denote the complex linear operator 

B = + ^ - Six)F'{{C, 0)) = [ (4.10) 
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where 



Di = -— + u-S{x)[a + b], 



L)2 = - 

The system (14. 8 p then reads as 



(4.11) 



D2 = -^ + uj-5{x)a. 



xix,t) = Cxix,t), C:=r'B= (4.12) 

Theorem 12.21 generahses to the equation (I4.12p : the equation admits unique solution x{x,t) G 
Cb(lR, H^) for every initial function x{x, 0) = xo ^ H^. Denote by e*^* the dynamical group of 
equation (14. 12p acting in the space H^. 



5 Laplace transform 

^ ^ POO 

Equation (I4.12p can be solved by the Laplace transform x{x,^) ■= / e^^^x{x^t)dt. The 

Jo 

Laplace transform is analytic function in the complex halfplane Re A > with the values in 
since the solution is bounded in H^. This implies that the resolvent R(A) := (C — A)~^ is also 
analytic for Re A > 0, with values in the space of bounded operators on H^. From the inversion 
of the Laplace transform we obtain 

ioo 

e*^* = / e^*R(A + £)rfA, t > 0, (5.1) 

2'ni J 

—ioo 

for any e > 0, where the integral converges in the sense of distributions of t G IR. 

We assume that the spectral condition (13.50 holds from now on. Then the resolvent admits 
analytic continuation from ReA > to the complex plain with the cuts C+ = [iuj,ioo), C_ = 
(—zoo, —iuj], and with the pole of order two at A = as detailed in Appendix A. Furthermore, 
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for A G C+ U C_, the resolvent R(A ± e) has right and left limits R(A ± 0) as £ — 0. Then (15.11) 
implies that 

e^* = — — [ e^*R(A) d\ - — f e^*(R(A + 0) - R(A - 0)) dX, for any r e (0, u) 

\X\=r c+uc_ 

(5.2) 

by the Cauchy theorem. Setting t = 0, we obtain that 

1 = — ^ / R(A)rfA-— / (R(A + 0) - R(A - 0)) rfA = P° + P", (5.3) 
2111 J 2ni J 

\X\=r C+UC- 

where P° and P'^ stands for the corresponding Riesz projectors (see [37]) onto, respectively, 
the generalised null space of C , and onto the continuous spectral subspace. We will show in 
the next section that P° is the symplectic projection, and therefore, P'^ is also the symplectic 
projection. The projectors P", P'^ commute with C and with the group e^*. Let us note that 

pOgCt ^ — L / e^*R(A) dX, P'^e^* = — - f e^*(R(A + 0) - R(A - 0)) dX. (5.4) 

\\\=T c+uc_ 

The first equation holds since both sides are one-parameter groups of operators , and their 
derivatives at t = coincide. The second equation follows from (15. 2p and the fact that 1 = 
po ^ pc Therefore, becomes 



6 Invariant subspace of discrete spectrum 

Here we prove that P*^ is the symplectic projection onto the tangent space of the solitary 
manifold S at the solitary wave e^^ip^i- The real form of the solitary wave is e^^^^j where 
= ii'ujix), 0). The tangent space to S at the point e-'^^^ with parameters u, 9 is the hnear 
span of the derivatives with respect to 6 and u cf. Remark [3. 3 p i.e. 

T^^S = linear spa.n^je^^^uj{x), e^^duj^uj{x)y 

Notice that the operator C corresponds to 6' = since we have extracted the phase factors e*^ 
from the solution in the process of linearization (14. ip . The tangent space to S at the point $^ 
with parameters (a;, 0) is spanned by the vectors 

To(a;):=j$^, T,{lu) := d^^^. (6.1) 

Observe that (12.71) and its derivative in u give the following identities: 

DaV'a. = Bi{d^^^) = (6.2) 

These formulae imply that the vectors Tq and Ti lie in the generalised null space of the non- 
self-adjoint operator C defined in (I4.12p and in fact Theorem IA.6l u) implies: 

Lemma 6.1 Let the spectral condition Ii3.5\) hold. Then the generalised null space of C is two 
dimensional, is spanned by Tq,Ti, and 

CTo = CTi = To. 
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We also introduce the symplectic form for the real vectors i/j and rj by the integral 

ai^,,) = IUM'i^ = l('P,,.-M, (6.3) 

where (•, •, ) stands for the scalar product in IR^. Since a' ^ a/C'^ then by Lemma [3^ 

= -^iTo,T,) = U^j \^^\^dx ^ 0. (6.4) 

Hence, the symplectic form VL is nondegenerate on the tangent space Tuj,o^i i-^- T^ qS is a 
symplectic subspace. Therefore, there exists a symplectic projection operator from L^(IR) onto 

Lemma 6.2 The operator Y*^ , defined in \5.^) . is precisely the symplectic projector from L^(IR) 
onto T^flS, and, furthermore, it may be represented by the formula 

P'^^p = boTo + biTi with — ^Jdq = ,Ti) , fi^^bi = il{^jj,To). (6.5) 

Proof The coincidence of both definition (15.31) and (16. 5p of operator P" follows by the Cauchy 
residue theorem from the formulas (17. 2p - (17.41) for the resolvent. ■. 

Corollary 6.3 P^ = 1 — P'' is also symplectic projector. 

Remark 6.4 Since To{uj),Ti{uj) lie in H^{]R) the operator P" extends uniquely to define a 
continuous linear map H~^(JR) — > T^^S , which is still designated P". In particular this operator 
can be applied to the Dirac measure 5{x). 

Using the Taylor expansion for the e'^* at A = and the identity AR(A) = CR(A) — 1, we 
obtain by (l5^ 

pOe*^* = (1 + Ct)P° (6.6) 



Remark 6.5 On the generalised null space itself C"^ = by Lemma \6. 1\ and so the semigroup 
gtc j-^dnQQg to 1 + Ct as usual for the exponential of the nilpotent part of an operator. 

7 Time decay in continuous spectrum 

From formulas (15. 5[ (16. 6p we see that the solutions x(^) = e*^*Xo of th^ linearized equation 
(14.12P do not decay as t — oo if P°xo 7^ 0. On the other hand, we do expect time decay of 
P'^x{t)j a consequence of the Laplace representation (15. 4p for P'^e*^*: 



pc^Ct 



— I e^*(R(A + 0) - R(A - 0)) c/A. (7.1) 

Tii J 



27fi 
c+uc 



The decay for the oscillatory integral is obtained from the analytic properties of R(A) for 
A G C+ U C_. The resolvent R(A) is an integral operator with matrix valued integral kernel 

R(A, X, y) = r(A, x, y) + P(A, x, y), (7.2) 
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where the columns of matrices F and P are given in (]A. 24p . (lA. 25|) . (lA. 271) . (lA. 280 : 

r(A,x,y) 





/ 


1 


1 


\ 








4A;_ 








i 


i 






\ 




4A;_ 


/ 


1 


/ gifc 






2D 




+ kl 


— 



jfc- |x| 



«a — 2k^ 

-tl3 



i(3 

-ia + 2A; I 



gjfc_|y| l^ik-\y\ 



(7.3) 



(7.4) 



Here fc±(A) = \J—uj =F ^A is the square root defined with cuts in the complex A plane so that 
A;±(A) is analytic on C\C± and ImA;-|-(A) > for A G C\C±. The constants a, /3 and D = D{\) 
are given by the formulas 

a = a + 6/2, (3 = 6/2, D = 2ia{k+ + A;_) - Ak+k^ + - (3^. 

Recall from Section |X] that D{\) 7^ for A G C+ UC_. Clearly in order to understand the decay 
of P'^e*'", it is crucial to study the behaviour of R(A,a;,?/) near the branch points A = ±20; 
(where k± vanish). 

We deduce time decay for the group P'^e**" by means of the following version of Lemma 10.2 
from [18], which is itself based on Zygmund's lemma [lH p. 45]. 

Let JF : [0, 00) ^ B be a function with values in a Banach space B. Let us define the 
B-valued function 



I{t) = le-''^T{0 dC 



Lemma 7.1 Suppose that J-'{0) = 0, and for some 6 > 

J^" e L\6, 00; B), (7.5) 

and 

r'(C) = OiC), CiO (7.6) 
in the norm of B for some p G (0, 1). Then I{t) G Cb{s, 00; B) for any e > 0, and 

I{t) = 0{t~^'P) as t^oo 

in the norm of B. 

For /5 > 2 let us introduce a Banach space which is the subset of distributions which are 
linear combinations of Ljj functions and multiples of the Dirac distribution at the origin with 
the norm: 

U + C5ix)\\M,:=mLl + \C\. (7.7) 

We will apply Lemma [TTT] to the function J-'{X) = R(A + 0) — R(A — 0) with values in the Banach 
space B = B^Aip, L'^p) , the space of continuous linear maps ^ L'^p for any [3 >2. 

Theorem 7.2 Assume that the spectral condition \3.5^) holds so that A = zs the only point 
in the discrete spectrum of the operator C = C{uj). Then for (3 >2 

\\P%^% = 0(^3/2^, t ^ 00. (7.8) 
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First we use the formulas (17.11) and (17. 2p to obtain 

-27riP^e*^*= J e^*(r(A + 0)-r(A-0))rfA + J e^\P{X + 0) - P{X - 0)) dX (7.9) 



c+uc_ 



c+uc_ 



Next we apply Lemma 17.11 to each summand in the RHS of (17. 9p separately. Then Theorem 
17.21 immediately follows from the two lemmas below. 



Lemma 7.3 If the assumption of Theorem \ 7. S\ hold then 

J e^\T{X + 0)-T{X-0)) dX = 0{t-^/^) 



t ^ oo 



(7.10) 



c+uC- 



in the norm B. 



Proof We consider only the integral over C+ since the integral over C_ can be handled in the 
same way. The point A = ia; is the branch point for k+, therefore, if A G C+ then since k- is 
continuous across C+ 

r(A + 0) - r(A - 0) = r+(A + o) - r+(A - o), 

where F"*" is the sum of those terms in F which involve /c+. Let us consider, for example, F^^. 
The expression (17. 3p implies for y > that 



rn(A,a;,?/) 



0, 



^ik+x/^-ik+y _ ^ik+y\ 



Ak. 



X < 0, 



< X < y, 



X > y. 



For A G C+, the root k+ = \J —uj — iX is real, and k+{X + 0) = — A;+(A — 0). Then, for y > 

sm{y^\x\) sm{^/C\y\) 



F+(A + 0,x,i/)-F+(A-0,x,2/) 



-e(x) 



(7.11) 



where ( = —uo — iA, and Q{x) = 1 for x > and zero otherwise. The second derivative of the 
fianction /(C) = ""^^'"'^r^^'^'^ satisfies 



ir(c)i 



sin(v^|a;|)sin(v^|y|)(|xJ^J^M^ ^ 2cos(v^|a;|) cos(v^|?/|)|a;||i/| 



4Cv^ 



siYi{^\x\) cos(v/C|y|)|y| + cos {^\x\) sin(v/C|y|)|a^l , ^ C(l + |xp)(l + |yp) 



2C 



For y < an identical calculation leads to the same bound. Therefore the operator valued 
function J^{() = Tii{X + 0) — Tii{X — 0) satisfies the conditions (17.50 and (17.60 of Lemma mi 
with ( = —uj — iX, p = 1/2 and B = i3. ■ 
Next we consider the second summand in the RHS of (17. 9p . 
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Lemma 7.4 In the situation of Theorem\7.2 



J e^\P{X + 0)-P{X-0)) dX = 0{t~^/^), (7.12) 
c+uC- 

in the norm B. 

Proof We consider only the integral over C+ and one component of the matrix P, for example, 
Pii: 

(ia- 2A;_)e**^+(l^l+ls'l) + {-ia + 2A;+)e^^-(l^l+IJ'l) + ip(e'^-\y\+ik+\x\ _ ^ik+\y\+ik.\x\^ 
Pii{X,x,y} - 2ia{k+ + k.) - 4k+k. + - " 

Denote ( = —u — iX, then k^ = ^/(, A;_ = a/— 2u; — (. A Taylor expansion in ^/( as C — ^ 
0, Im(^ > implies 

Pn(A, X, y) = Po + Pi{x, y)^'^ + P^i^x, y)0{C). 

where \Pj[x,y)\ < Cj{l + |a;p)(l + |yp), j = 1,2. Therefore, if A G C+ then 

J'iC) = ^ii(A + 0) - Pn(A - 0) = 0(C^/'), C - 

in the norm of B. Similarly, differentiating two times the function Pii{X,x,y) in A, we obtain 
that 

-^"(C) = -A"i(A + 0) + Pri(A-0) = O(r=^/2)^ ^,J = l,2, c^o 

in the norm of B. Moreover, J-'"{C) ~ C~^^^ as C ^ oo- Therefore, the function J^(C) satisfies 
the conditions (17. 5p and (17.61) of Lemma mi with p = 1/2 and B = i3. ■ 

8 Bounds for small t 

As a starting point for the method of majorants in Section [10] we will need also some estimates 
on the group e*^* for small t. First note that the function e^*Xo belongs to Cb(]R, H^). This 
follows from a theorem analogous to Theorem 12.21 for solutions e*^*Xo of the linearized equation 
(I4.12p . with initial condition xo ^ H^. Moreover, energy and charge conservation imply that 

We^'xoWm <c\\xo\\m, teM. (8.1) 

For a further application in section 111.41 we need a bound for the action of e^^ on the Dirac 
distribution 6 = 6{x). Thus let xs{x,t) be the solution to the linearized equation (14.41) with 
X5{x,0) = 6{x) and e^^6 its real vector version. Note that, by Theorem 17. 2[ we have e^*(5 G 
e, oo; L°^p), for every e > 0, and (3 > 2. The next lemma gives the small t behaviour: 

Lemma 8.1 The following bound holds 

||eC*5||ioo = C(ri/2), t^O. (8.2) 
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Proof By the Duhamel representation for the solution to (14.41) . we obtain 

t 

X8{x, t) = W^{t)5 - / ds(ax5(0, s) + 6Re (^^(O, s))w^{t - s)5 (8.3) 



where a and h are defined by (14. 9p . and W^jif) is the dynamical group of the modified Schrodinger 
equation 

ix{x, t) = -x"ix, t) + uxix, t). (8.4) 

Note that 

W^{t)5 = -jL= e*^-*^* (8.5) 

Therefore (18.31) becomes 

t 

1 2 r ^ . -2 



X^(x,t) = -==e%"-^- / e^^--"(*--^(ax^(0,g) + 6Re(x^(0,g))cig (8.6) 

V47r«t ^ J47ri(t-s) ^ ^ 



Denote '^{x,t) = \ft X5{x,t). Then 

t 

<;{x, t) = ^== e'^-'"* - / , ^ e'^-'"^*-^^ (a?(0, s) + 6Re (?(0, s))ds. (8.7) 



^47r? ^ ^47ri(t - s)s 

Therefore, 

11^ } 1 

IkWIU- < + 7;v^(|a| + 1^1) / ; |k(s)||iooC?s, t > 0. 

^\/7r / J or. (-1- — q\q 







2\/7r 2 y 7r^(t-s)s 

7rw'(t- s)s 



(is 

Since / — , = 1, we obtain the bound 



1 1 

k(t) koo < — — --= 

- 2v^l-iv^(|a| + |6|) 
if t is sufficiently small. ■ 

9 Modulation equations 

In this section we present the modulation equations which allow a construction of solutions 
ip{x,t) of the equation (11.11) close at each time t to a soliton i.e. to one of the functions 

Ce^^-v^l^l, C = C{uj)>0 

in the set S described in section |3] with time varying ("modulating") parameters {uj,6) = 
{uj(t),9(t)). It will be assumed that ip{x,t) is a given weak solution of (11.11) as provided by 
Theorem 12.21 so that the map t —>■ ,t) is continuous into H^{]R). The modulation equations 
follow from the ansatz for the solution which is explained next. Recall that we defined 

$,(x) = (Ce-v^l^l,0) = (^^,0) (9.1) 
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so that '4'{x,t) = e^^^*^^uj{t)ix) is a solution of (14. 7p if and only ii 9 = cu and ci; = 0. Here it is 
to be understood that C = C{uj{t)) is determined from ujit) via (13.31) . We look for a solution 
to (14.71) in the form 

ijix,t) = e^'^'\^^^t)ix) + xix,t)) = e^'^'^^{x,t), ^(x, t) = <l>^(i)(x) + t). (9.2) 

Since this is a solution of (14.71) as long as x = and 6 = ui and ci; = it is natural to look for 
solutions in which x is small and 

e{t) = f uj{s)ds + -iit) 
Jo 

with 7 treated perturbatively. Observe that so far this representation is underdetermined since 
for any {uj{t),6{t)) it just amounts to a definition of x] it is made unique by restricting x{t) 
to lie in the image of the projection operator onto the continuous spectrum = P'^(u;(t)) or 
equivalently that 

P°x(t) = 0, P? = P°(^(t)) = / - P'^iuit)) (9.3) 

(The projection operators are as defined in (15. 3p and ( 16. 5p ). An equivalent formulation of (19. 3p 
is to say that e^^x is required to lie in the symplectic normal space N^(t),e(t)S. This is equivalent 
to imposition of the following orthogonality conditions (at each time t): 

fi(x(t),To(a;(t)) = = fi(x(t),Ti(a;(t)), (9.4) 

where fl is the symplectic form introduced previously. Writing x{t) = (xi(^)) X2it)) the orthog- 
onality conditions reduce to 

J Xi(x,t)Ce"v^l^lda; = 0, and J X2(x, t)5^(Ce-^l^l) = 0. (9.5) 

Now we give a system of modulation equations for uiit), •yit) which ensure the conditions 
(19. 5p are preserved by the time evolution. 

Lemma 9.1 (i) Assume given a solution of d^^. 7| ) with regularity as described in theorem \2.2\, 
which can be written in the form Ii9.^) - I{9.3\) with continuously differentiable uj{t), 9(t). Then 

X = Cx- djd^^^ + 7J-^($. + X) + Q (9.6) 

where Q(x,u;) = -6{x)r\F{<l>^ + x) -F{<l>^) ~F'{<^^)x), and 

(9A-9^P0X,^) ^ ' 

. ^ (jP°(g^i>^-g^P°x),P°Q) . . 

where P° = P°(u;(t)) is the projection operator defined in ^6. 5]) and d^iP^ = 9^P°(u;) evaluated 
at u = uj{t). 

(a) Conversely given ip a solution of ^-ID as in theorem \2.2\ and continuously differentiable 
functions uj(t), 6(t) which satisfy \9. 7| )-( fgr^) then x defined by ^9.2\) satisfies ^9.d^) and the 
condition holds at all times if it holds initially. 



14 



Proof This can be proved as in [6, Prop. 2. 2]. ■ 
It remains to show, for appropriate initial data close to a soliton, that there exist solutions 
to (I9.7l) - (l9.8p . at least locally. To achieve this observe that if the spectral condition (13.51) holds 
then by Lemma 13.41 the denominator appearing on the right hand side of fl9.7p and (19.81) does 
not vanish for small This is because 

{d^^s.^s) = \d^j \^.?dx ^ (9.9) 

as discussed in section 121 This has the consequence that the orthogonality conditions really 
can be satisfied for small x because they are equivalent to a locally well posed set of ordinary 
differential equations for t {6{t),uj{t)). This implies the following corollary: 

Corollary 9.2 (i) In the situation of (i) in the previous lemma assume that h9. 9\) holds. If 
WxWl^p '^s sufficiently small for some p, (3 the right hand sides of ^9. 7| ) and l{9.8\) are smooth in 
9,u and there exists continuous IZ = 7l{uj,x) such that 

|7(t)| <7^|x(o,t)|^ \djit)\<n\xio,t)\'. 

(a) Assume given ip, a solution of ( [y^. 7[ ) as in Theorem \2.S[ IfujQ satisfies i\9. 9\) and x(x, 0) = 
g-jfo^(3;^Q) _ <|)^^(x) is small in some norm and satisfies \9.'J\) there is a time interval on 
which there exist functions t ^ {uj(t),'y(t)) which satisfy ^9. 1\ )- f97^) . 



10 Time decay for the transversal dynamics 

In this section we state our Theorem 110.21 on the time decay of the transversal component x{t) 
in the nonlinear setting, leaving the proof to the next section. Theorem 110.21 will be used to 
prove the main theorem in Section [121 First we represent the initial data ipo in a convenient 
form for application of the modulation equations: the next Lemma will allow us to assume that 
(19.31) holds initially without loss of generality. 

Lemma 10.1 In the situation of Theorem 3.1 there exists a solitary wave ijjcoQ = Coe"^^!^! 
satisfying the spectral condition 1^3. 5\) such that in vector form 

M^) = e^'^"($^o(x) + Xo{x)), = (V^^o, 0), 

and for Xoi^) we have 

P°(cuo)(xo)=0, (10.1) 

and 

WxoWLlnm = d = 0{d) asd-^0. 
Proof By (19.41) . the condition (110. ip is equivalent to the pair of equations 

n{e-^''iJo - $^„,To(cDo)) = 0, nie'^'^^o - <^>^o> ^i('^o)) = 0, 

where Tq^lj) = j^uj, Ti{uj) = d^^^^. For ipo sufficiently close (in L^) to e^^°^ujo the existence 
of '^o follows by a standard application of the implicit function theorem if we show that the 
Jacobian matrix 



:i0.2) 
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with ipQ = e^^°^i^Q, is non-degenerate at u = ujq and 9 = 6q. But this is equivalent to the 
non- degeneracy of the matrix 







;io.3) 



which holds by (EJl). ■ 
In Section [12] we will show that our main Theorem 13.71 can be derived from the following 
time decay of the transversal component xi^)'- 



Theorem 10.2 Let all the assumptions of Theorem 3.1 hold. For d sufficiently small there 
exist functions t ^ {uj{t),'y{t)) defined for t > such that the solution ip{x,t) of "A ) ca?^ 
be written as in l{9.^9.3\) with ^9. 1^9. 8\) satisfied, and there exists a number M > 0, depending 
only on the initial data, such that 

M{T)= sup + +(i + t)3(|^| + |^|)|] (10.4) 

0<t<T 

uniformly in T > 0, and M = 0{d) as d ^ 0. 

Remarks 10.3 (0) This theorem will be deduced from Proposition \ll.l\ in the next section. 

(i) Theorem \2.^ implies that the norms in the definition of M are continuous functions of 
time (and so M is also). 

(a) The result holds also for negative time with appropriate changes since ilj{x,t) solves 
( li.ij) if and only if ip{x, —t) does. 

(Hi) The result implies in particular that t^\6 — uj\ + t^\uj\ < C , hence uj{t) and 9{t) — tuj^ 
should converge as t oo while ~ e-'^^*''$a;(t)(a;) have limit zero in L'^piWl). 

(iv) The notation x{t) indicates the function x ^ xi^^t) ^■^ usual. 

11 Proof of transversal decay 

11.1 Inductive argument (proof of Theorem 110^2]) 

Let us write the initial data in the form 

Mx)=e'''i^^oi^)+Xoix)). (11.1) 

with d = llxollLinHi sufficiently small. By Lemma [10.11 we can assume that P°(ci;o)(xo) = 
without loss of generality. Then the local existence asserted in Corollary 19.21 implies the 
existence of an interval [0,ti] on which are defined functions t i— > {uj{t),'y{t)) satisfying 
(I9.7I) - (I9.8I) and such that M(ti) = p for some ti > and p > 0. By continuity we can make p 
as small as we like by making d and ti small. The following Proposition is proved in section 
[m below. 

Proposition 11.1 In the situation of Theorem \10.2\ let M{ti) < p for some ti > and p > 0. 

Then there exist numbers di and pi, independent of ti, such that 



M(ti) < p/2 (11.2) 



ifd= IIXo||Lin/fi < di and p < pi. 
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Assuming the truth of Proposition 111.11 for now Theorem 110.21 will follow from the next argu- 
ment: 

Consider the set T of ti > such that {uj(t),'y(t)) are defined on [0, ti] and M(ti) < p. This set 
is relatively closed by continuity. On the other hand, flll.2p and Corollary 19.21 with sufficiently 
small p and d imply that this set is also relatively open, and hence supT = +oo, completing 
the proof of Theorem I1U.2I ■ 

11.2 Frozen linearized equation 

A crucial part of the proof of Proposition 111.11 is the estimation of the first term in M, for 
which purpose it is necessary to make use of the dispersive properties obtained in sections E] 
and [71 Rather than study directly (19. 6p . whose linear part is non-autonomous, it is convenient 
(following [31 15) to introduce a second ansatz, a small modification of (19. 2p . which leads to an 
autonomous linearized equation. This new ansatz for the solution is 

ij{x,t) = e^^^M + e~^^^-~^\), where e{t) = uit + ^o, = andcji = uiti) (11.3) 

so that, T] = e^^^~^'^x ^-nd x = e~^^^~^\. Since 

equation (19. 6p implies 

f, = j-\cui - Lu)r] + e'^'''^C{e-^^'-'\) + e'^'-'^ {j'^i^u. - ^d^^. + Cl[e-'^'-'^i\). (11.4) 
The matrices C and e-''^, where (j) = 6 — 0, do not commute hence we need the following lemma: 



Lemma 11.2 



Ce^^ - e^^C = 6{x)bsm(f) a, where = ( J ^ ^ = 2a'Cl (11.5) 



Proof 



Ce^'f' - e^'t'C = i D2 W COS0 -sin0 \ / COS0 -sin0 \ / D2 \ 

Y — Di y \^ sin0 cos0 J y sin0 cos0 J y — Di J 

_ / (D2 -Di)sin0 \ _ f 5{x)b sin (p ^ \ m 

~\ (Di -D2)sin0 J ~ -(5(x)6sin j ' 

Using Lemma [11.21 we rewrite equation (111.40 as 

?) = j-\uji - uj)r] + Cr] + e^'^^-'^H -6{x)b sm{e - e)ar, + r^^^^ - Cod^^^ + q,[e-''^^-'^\]). 



To obtain a perturbed autonomous equation we rewrite the first two terms on the RHS by 
freezing the coefficients at t = ti. Note that 
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where V = a + bPi, Vi = V{ti), and Ci = C{ti). The equation for r] now reads 

7) = Cir]-j-^6{x){V-Vi)r]+e^'^^'^^(^-6{x)bsm{e-e)a7]+r^^-iod^^^+Q[e-^^^-^^ (11.6) 

The first term is now independent of t; the idea is that if there is sufficiently rapid convergence 
of uj{t) a.s t —>■ oo the other remaining terms are small uniformly with respect to ti. Finally the 
equation (111.61) can be written in the following frozen form 

ri = Cir] + h (11.7) 

where 

fi = -j-^6{x){V~Vi)r] + e^^^~^^(^~6{x)bsm{e-e)a7] + j-^<l>-djd^<^^ + Q[e-^'^^-^^ (11.8) 



Remark 11.3 The advantage of Ijll-T] ) over Ii9.6\) is that it can be treated as a perturbed 
autonomous linear equation, so that the estimates from section can be used directly. The 
additional terms in fi can be estimated as small uniformly in ti: see lemma 11. 4 below. This 
is the reason for introduction of the second ansatz lill.3\) . 

Lemma 11.4 In the situation of Proposition [77771 there exists c> 0, independent ofti, such 
that forO<t<ti 

\a{t) - ail + \b{t) - fell + \9{t) - e{t)\ < cp. 

Proof By ([I02D 

sup il + t')i\^it)\ + \djm<Mit,) = p. (11.9) 

0<t<ti 

Therefore 

\a{t) - a{ti)\ = I / d(r)rfr| < c( sup (1 + T^)\dj{T)\) f < cp, 

since |a(T)| < c|u;(r)|. The difference \b(t) — b(ti) \ can be estimated similarly. Next 

e{t)-e{t) = /* cu(r)rfr + 7(t) -cu(ti)t- 7(0) = f\uj{T) - uj{ti))dT + f' i{r)d(fil.lO) 
Jo JO JO 

ft ftl ft 
= — / uj{s)dsdT + / 7(r)(ir. 

JO Jt jo 

By (111.91) the first summand in RHS of flll.lOp can be estimated as 

ft ftl 1 

<c sup (1 + s)^+^|c:;(s)| / / ——dsdT<cp 



0<sW ' ' ' "Jo Jr (1 + S)2+^ 

since the last integral is bounded for t G [0,ti]. Finally, for the second summand on the RHS 
of flll.lOp inequahty (fTOj) implies 

d,T 

< cp 



'1 T 

ft f ^J- 

/ 7(r)c/r| <c sup (l + r')|7(r)| / - 

-'0 0<T<tl J i 



+ r2 

t 
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11.3 Projection onto discrete and continuous spectral spaces 

From sections [6] and [7| we have information concerning U{t) = e^^*, in particular decay on the 
subspace orthogonal to the (two dimensional) generalized null space. It is therefore necessary 
to introduce a further decomposition to take advantage of this. Recall, by comparing fl9.2p and 
firr^ that 

^ = e^i^^-^^X and P°x(t) = (11.11) 

Introduce the symplectic projections = and PI = Pj onto the discrete and continuous 
spectral subspaces defined by the operator Ci and write, at each time t G [0,ti]: 

V{t) = g{t) + h{t) (11.12) 

with g{t) = P5?7(t) and h{t) = P5?7(t). The following lemma shows that it is only necessary to 
estimate h(t). 

Lemma 11.5 In the situation of Proposition [7777]. assume 

sup (|cu(t)-cui| + |0(t)-^i(t)|) = A 

0<t<ti 

is sufficiently small. Then for < t < ti there exists a linear operator S(t), bounded on 
L'^p n , and c(A,co'i) such that ri{t) = 'E{t)h{t), and 

c{A,uJiy^\\\h\\Loo^r,m < ||^?||L-^n//i < c{A, uJi)\\\h\\Loo^nHi ■ (11.13) 

Proof Explicitly we write 

Vit) = hit) + git), git) = 6o(t)To(^i) + 6i(t)Ti(a;i) (11.14) 

where 6o)^i are chosen at each time t to ensure that f2(/;,(t), ro(ct;i)) = f2(/;,(t), ri(ct;i)) = 0. 
Using the fact that (since P°x(t) = 0) 

nie-^^'~'\Muoit)) = = nie~'^'-'\,Tiiuoit)) 

this means that h^, hi are determined by 

-li^Mt) =niv{t),Tiiui)) = Qi7]it),TiiuJi)-e^^'-'^Tiiujit))) (11.15) 
^iu:M{t) =fi(r/(t),To(cui)) =^](r/(t),To(cUl)-e^■(^-^')To(cu(^))). (11.16) 

From these it follows that there exists c > such that ||5'(t) nJ/i < cA||?7(t)||ioo j-^j^i and 
hence (I11.13P follows as claimed. ■ 

11.4 Proof of Proposition 111.11 

To prove Proposition 1 1 1 . 1 1 we explain how to estimate both terms in M, (110. 4p . to be < p/4, 
uniformly in ti. 

Estimation of the second term in M. As in Corollary 19.21 we have 

MitV 

m)\ + Ht)\ < co|x(0,t)|' < co-^^-^, t < ti, 
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since |x(0,i(:)| < ||x(t)||L°o^. Finally let pi < l/(4co) to complete the estimate for the second 
term in M as < p/4. 

Estimation of the first term in M. By Lemma [11.51 it is enough to estimate h. Let us apply 
the projection to both sides of (111.71) . Then the equation for h reads 

/i = Ci/i + P^fi (11-17) 

Now to estimate h we use the Duhamel representation: 

h(t) = U(t)h(0)+ f'u(t-s)PlU{s)ds, t<ti. (11.18) 
Jo 

with U{t) = e*^^* the one parameter group just introduced. Recall that P?/i(t) = for t G [0, ti]. 
Therefore 

\\U{t)h{0)U^^<c{l + t)-'/^h{0)^,,^.< (11.19) 

by Theorem 17.21 and inequalities (18.11) and (111.131) . Let us estimate the integrand on the right- 
hand side of (|TLT8ll . 

Lemma 11.6 The integrand in U1.18\) satisfies the following bound: 

||[/(t-s)P^fi(s)|U^^ < c ^^ _ ^y/^li+t^ ,) {Ms)\\1^+pMs)\\l^^^, 0<s<t. (11.20) 

Proof We consider two different cases : t — s > u, and < t — s < u, where u = z/(a, b) is 
defined in Lemma 18.11 

i) t — s > u : We use the representation (I11.8P for fi and apply Theorem 17.2^ Corollary 19.21 
and Lemma 111.41 to obtain that 

\\U{t-s)Plh\\L^^ < c{u){l+t-s)-'/'\\Pl{i,m\M 

< ciu)ii+t-sr'^'(\vio,t)\'+pmt)\^ 

< c(z/)(l+t-s)-=^/2(^||r;(t)||i^^+ph(t)iU^^), t<ti. (11.21) 

ii) < t — s < z/ : We denote Q = 6{x)Q, and represent fi(x, s) as 

fi(x, s) = pis)5ix) + qis)^^ + r{s)d^^^. (11.22) 

where 

p{s) = -3-\V - Vi)r]{0, s) + e^^^-^^{bsm{e - e)cxr]{0, s) + Q[e-^(^-^')r/(0, s)]), 
is an IR^ valued function of time, and 

q{s) = -e-^^'~'^n, r{s) = -e-'^'~'^u 

are (2 x 2) matrix valued functions of time. Writing || ■ || for both the standard Euclidean and 
operator norms on these, we have, by Lemma lll.4[ 

|b(^)|| < c(\v{0,s)\' + p\v{0,s)\) < c(||r^(.)||i^^ +p||r/(.)|U^^) 
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and by Corollary 19.21 

||g(s)||,||r(.)||<c|r/(0,.)p<c||r/(.)||i.. 
Applying projector PI to fi we obtain 

P^fi(a;, s) = pis)6ix) + ^(s)*^ + r{s)d^^^ - P%{x, s). (11.23) 

By Lemma 18.11 for sufficiently small u we obtain 

\\U{t - s)p{s)6{x)\\,^^ < \\U{t - s)p{s)S{x)h^ < c{u)\\p{smt - 

< c(i.)(t-s)-^/'(h(s)||ioo^ + ph(s)|U«=^), 0<t-s<u. (11.24) 

By inequality (18.11) we have 

II Uit - s) {q{s)^^ + r(s)9^$^) < c||t/(t - s) [q{s)^^ + r{s)d^^^) \\hi (11.25) 

< c(||g(3)||||<l>J|Hi + (||r(s)||||aA||//0 <c||r7(s)||i^^, 0<t-s<i.. 

The definition (16. 5p of the projector implies immediately that 

l|P?filki<c(|b(.)|| + ||g(.)|| + ||r(.) 
Then, similarly to (111.251) . we obtain 

\\Uit - 5)P?fi|U^, < ciUs)\\l^^ + pMs)\\l^^), Q<t-s<u (11.26) 
Finally, ([IL23D-([IL26D imply 

\\U{t - 5)P?fi||L-^ < c(t - s)-"\Ms)\\l^^,^ + p||r/(s)|U^J, Q<t-s<v. (11.27) 



From ffTT:2TD and ffTTTTI) inequality ^TLW) follows. 
Now ffTTTS]) . (HHHD, ffTLTQl) and ffTT:20D imply 



t 

\\m\L-^<c{l+t)-'''\\m\\Llr.m+c, I _ ^y/,^^ ^ ^ _ (h(^)lli^^ + pMs)\\l-^) 
Multiply by (1 + t)^/^ deduce 



{l + tm\vmL-^<cd + c, J (^^^.^^(iVV- + .)1r^(.)||i^/. (11.28) 

since ||??(0)||/,inj:^i < d. Introduce the majorant 

m(t) := sup(l + s)^/^7]{s)\\l^ t < 

m " 
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and hence 

j fl _L +^3/2/1 I N-3 j (I ^f]3/2(i , \-3/2 

,n(t) < ai + c.,n'(0 / ^ ^ ^.^^ ds + / |^ ^ ^^^^ ^ f_ ,s. (11,29) 

It easy to see (by splitting up the integrals into s < t/2 and s > t/2) that both these integrals 
are bounded independent of t. Thus (111.291) implies that there exist c, C2, C3, independent of ti, 
such that 

m(t) < cd + pc2m{t) + C3m'^{t), t < ti. 

Recall that m{ti) < p < pi by assumption. Therefore this inequality implies that m{t) is 
bounded for t < ti, and moreover, 

m(t) < c^d, t <ti 

if d and p are sufficiently small. The constant C4 does not depend on ti. We choose d in (13.61) 
small enough that d < p/^Ac^). Therefore, 

sup(l + tf/lry(t)|U^ <p/4 

[o,ti] " 

if d and p are sufficiently small. This bounds the first term as < p/4 by (111. lip and hence 
M{ti) < p/2, completing the proof of Proposition lll.il ■ 



12 Soliton asymptotics 

Here we prove our main Theorem 13.71 using the bounds (110. 4p . For the solution ilj{x,t) to (11.11) 
let us define the accompanying soliton as s(x, t) = '(/'^(t)(a;)e*^W, where 9{t) = uj{t)+^{t). Then 
for the difference z{x,t) = il){x,t) — s{x,t) we obtain easily from equations (11.11) and (12.71) 

iz{x,t) = -z"{x,t) +^s{x,t) -iujd^s{x,t) - 5(x)(F(^(a;,t)) - F(s(a;,t))). (12.1) 

Then 

t 

z{t) =Wit)z{0) + I Wit ~ T)[ys{T) - tud^sir) - 5{x)[f{^{0,t)) - F{siO,T)))]dT, (12.2) 


where z(t) = z{-,t), s(t) = s(-,t), and W(t) is the dynamical group of the free Schrodinger 
equation. Since ■y(t) — 7+, uj(t) — uj+ = 0(t~^), and therefore 9(t) — u+t — 7+ = 0{t^^) for 
t — i> 00, to establish the asymptotic behaviour (13.71) it suffices to prove that 

z{t) =W{t)<l>+ + r+{t) (12.3) 

with some $+ e a(IR) n L\M) and \\r+{t)\\c,in)nL^u) = 0{t-^/^). Denote git) = 7s(t) - 
iujd^s{t), hit) = F{ij{0,t)) - F(s(0,t)) and rewrite equation (fTO) as 

t t 
z{t) = W{t)z{0) + W{t) f W{-T)g{T)dT - W{t) I W{-T)5{x)h{T)dT. (12.4) 



22 



Let us recall that uj{t), 'j{t) ^ t ^ as t ^ oo. Hence, for the second summand in RHS of (112.41) 
we have 

t oo oo 

W{t) J W{-T)g{T)dT = W{t) J W{-T)g{T)dT ~W{t) J W{-T)g{T)dT = W{t)(f)i + ri{t), 

t 

(12.5) 

where, from the unitarity in of the group W{t) and the t ^ decay of uj and 7, we infer that 

00 

(j)i = J W{-T)g{r)dT G H\ and ||ri(t)||^^i = ^(r^), t 00. 


gis2/(4i) 

Consider now the last summand in the RHS of (112.41) . Note that W{t)6{x) = — ^==, and 
\h{t)\ < c|x(0,t)| < c(l + t)-3/2 Therefore 

* 7 -ix2/(4r) y ixy{4{t-r)) 

W{t) W{-T)5{x)h{T)dT = W{t) / , h{T)dT- / , h{T)dT = W {t)(p2 + r2{t) . 

i i V-47rzr { JAm{t-T) 



(12.6) 

°?g-ixV(4r) 

Evidently, 02 = / —i= h{r)dT E Cb, and ||r2(t)||cb = 0{t'^), t 00. 

J \/—A'KiT 

^ 

Moreover, 02 ^ L"^, and ||r2(t)||L2 = C(t^^/^), t — 00. To see that this is indeed true change 
variable to r = 1/m in the definition to get: 

02(x) = , ^ r e"™"'/^ r]{u) du, 7]{u) = h{l/u)/u^'\ (12.7) 



a/— 47ri ^0 

Now hit) is bounded and it follows from the decay of h{t) that rj{u) is bounded as m — > 0. 
Therefore ri{u) is square integrable and so by the Parseval-Plancherel theorem 02 is square 
integrable as a function of y = x'^, and hence also as a function of x (since dy = 2xdx and 02 
is a bounded continuous function). Next we have r2{t) = —W{t)R{t) with 

R{x,t) = -jl= r^*e— V4 ^(^) = _A^F^^,.„an)ri{n), 
^f—Am Jo \J —A-n% 

where C,t{u) is the characteristic function of the interval (0, 1/t). The function r]{u) is bounded, 
hence HCt^^lU^ = ct~^/^ and therefore ||r2(t)||i2 = 0{t~^l'^), t 00. To conclude, using (112. 4p . 
ffT23|l . and we obtain ffTO) with 0+ = ;z(0) + 0i + 02 and r+{t) = ri(t) + r2{t). The 

t — > —00 case is handled in an identical way. ■ 



Remark 12.1 The expression ( 12. 1 ) for 02 as a Fourier transform implies immediately that 
|02|, and hence |$+| also, tend to as \x\ —>■ 00 by the Riemann-Lebesgue lemma. This same 
expression could be used with Zygmund's lemma to obtain more detailed decay properties of 
02 and hence of . The decay rate would be determined essentially by the regularity of the 
function h{t) in addition to the decay rate of the initial data. 



A The resolvent 



A.l Calculation of the matrix kernel 

The derivation of the time decay of the solution to the linearized equation (14.41) in section [6] re- 
quired an analysis of the smoothness and singularities of the resolvent R(A) and its asymptotics 
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for A — > oo. Here we will construct its matrix integral kernel explicitly 



R(A,x,?/) 



which is the solution to the equation 



R2iiX,x,y) i?22(A,x,?/) 



iC-\)Ri\,x,y) = 5ix-y) 



1 
1 



(A. 1) 



(A. 2) 



Calculation of first column For the first column Rt(\, x, y) := f ^^^1^' ^\ I of the matrix 

' '^^ \^ R2i{X,x,y) J 

R(A,x, I/) we obtain 



iC-\)Rji\,x,y)=6ix-y) 







(A. 3) 



If a; 7^ and X 7^ (jSZH takes the form (cf. fliATl) . fl412D ) 

Ri{X,x,y) = 



-A 
-Di 



D2 

-A 



-A 



dx"^ 



\ 



LU 



-A 



Ri{X,x,y) = 0, x^O, x^y. 

(A. 4) 

The general solution is a linear combination of exponential solutions of type e^^'^v. Substituting 
into (lA. 41) . we get 

For nonzero vectors v, the determinant of the matrix vanishes, 



0. 



Then k'^ + u = =FiA. Finally, we obtain four roots ±A;-|-(A) with 



k±{X) = \l —uj =F iX, 



(A. 6) 



(A. 7) 



where the square root is defined as an analytic continuation from a neighborhood of the zero 
point A = taking the positive value of \m.^ —uj at A = 0. We choose the cuts in the complex 
plane A from the branching points to infinity: the cut C+ := [iu;,ioo) for /c+(A) and the cut 
C_ := [— icij, — ioo) for k_{X). Then 



Imfc±(A)>0, AgC\C±. 
It remains to derive the vector v = (f i, V2) which is solution to ( lA. 51) : 



(A. 



V2 



k\^ uj ±.iX 

■^^"^ = -X 



v\ = ±ivi. 



Therefore, we have two corresponding vectors ^± ~ |^ j- ■ J ^^'^ linearly independent 

exponential solutions 



V-e 



±ik—x 
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Now we can solve the equation (]A. 3|) . First we rewrite it using the representations fl4.12p and 
( 14. lip for the operator C, 



-A 



— u 



^21 (A, X,?/) 



-A 



5{x-y)i^\+5{x)i ° aUi?n(A,0,,) 



(A. 9) 

Let us consider y > for the concreteness. Then the RHS vanishes in the open intervals 
(— cxD, 0),(0, y) and (y, cxd). Hence, for the parameter A outside the cuts C±, the solution admits 
the representation 



Ri{\,x,y) 



B-e 



-ik+x„ 



B-e 



•ik-x„ 



-ik-X„ 



, X <0, 

+ B+e'''+'^v+ + file^^-^t;_, 0<x <y, (A. 10) 



X > y 



since by f lA. 81) . the exponent e~*^*^ decays for x — oo, and similarly, e*'^='='^ decays for 
X ^ oo. Next we need eight equations to calculate the eight constants . . . , C_. We have 
two continuity equations and two jump conditions for the derivatives at the points a; = and 
X = y. These four vector equations give just eight scalar equations for the calculation. 
Continuity at x = y: Ri{y — 0,y) = Ri{y + 0,y), i.e. 

BZv+/e+ + 5rf_/e_ + 5+i;+e+ + Btv^e_ = C+v+e+ + C_w_e_, 

where e± := e*'^*^. It is equivalent to 

B-/e+ + 5+e+ = C+e+, 

Bz/e- + Bte- = C_e_. 
Continuity at x = 0: Ri{—0,y) = Ri{+0,y), i.e. 

+ A_v_ = B-V+ + BZv_+ Blv+ + B^v. 

that is equivalent to 



(A. 11) 



f A+ = B-+Bt 
A_ = B- 



B+. 



(A. 12) 







Jump at X = y: R'j{y + 0,y) = R'i{y — 0,y) + \ ^ ) , where prime denotes the derivative in 



X. Substituting (lA. lOp . we get 



tA;+-B_i_t'+/e+ — ik^B_v^/e- + iA;+-B^t'+e+ + ik^B^V-C- 





-1 



(A. 13) 
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Noting that 



we get 





-1 



f + — V- 



ik_C_e_ = -ik_Bz/e_ + ik_B^e_ - 



(A. 14) 



(A. 15) 



After substituting of C± from ( lA. lip , the constants B^ cancel and we get 

e+ . e_ 



b: 



4k. 



Jump at X = 0: R'j{+0,y) = R'j{-0,y) - 



BZ 

a + b 
a 



Ak. 



(A. 16) 



Ri{-0,y). Substituting flA. lOp . we 



get 

- ik+BZ^v^ - ik_BZv^ + ik^B^v^ + ik_Btv^ = -ik+A+v+ - ik^A^v^ - M{A+v+ + A^v^), 

(A. 17) 



where M is the matrix 



a + b 
a 



. Note that 



^ a b 

where a = a -\ — , p = -■ 
2 2 



Mv- = av- + 
Then (lA. 17p becomes 

-ik^B~ + ik+B+ = -ik+A+ - A+a - A_f3, 

-ik^Bz + ik_B+ = -ik^A_ - A+(3 - A_a. 
Substituting here (1A. 12^ . we get after cancellations, 

' {2ik+ + q)BX + PBl = -aB- - PBz 

< 

^ (3B+ + {2ik_ + a)Bt = -/SB' - aBz 
Hence, the solution is given by 



(A. 18) 



B: 



1_ / 2ik_ + a 
D [ 



a P\i B^- 



—P 2ik^ + a I \ P a J 
where D is the determinant 

D := {2ik+ + a){2ik_ + a) - p^, 
and B^ , Bz are given by ( lA. 16p . The formulas (lA. 16p and ( lA. 19p imply 



(A. 19) 



(A. 20) 



2D 



2ik_a + 0^ - P' 
2kl 



e+ + z/3e_ , Bl 



2D 



' ^ 2ik+a + a^-p^ 
— «pe+ H — e„ 



2k. 



(A. 21) 
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Using the identities 

2ik_a + - (3'^ = D - 2ik+a + Ak+k-, 2ik+a + - (3'^ = D - 2ik_a + Ak+k. 
let us rewrite flA. 2ip as 



+ ^((^« - 2A:_)e+ + i/3e.), = ^ - ^(^/9e+ + - 2A;+)e_). (A. 22 



Ak^ 2D 



Finally, the formulas flA. lOD - flA. 12p . flA. 161) and flA. 221) give the first column Ri{X,x,y) of 
the resolvent for y > 0: 

RiiX,x,y) = Tj{X,x,y) + PiiX,x,y), (A. 23) 



where 



Ti{X,x,y) 



Ak 



1 '^ik+\x-y\ _ ^ik+i\x\ + \y\)^^^ _ _J_^^ik.\x-y\ _ ^ik-(\x\ + \y\)^^^ 



Ak. 



and 



Pi{X,x,y) = ^[((?a-2A;_)e'^+(l"l+l^l) + 2/?e'('=+l"l+^-l^l))t;+ 
- (z/5e'(fc-l-l+fc+l?/l) + (^a - 2A;+)e*^-(l"l+IJ'l))t;_ 



(A. 24) 



(A. 25) 



Calculation of second column The second column is given by similar formulas with the 
vector ^ Q ^ instead of ^ ^ ^ in flA. 9|) . Then ^ i ^ flA. 131) is changed by ^ q ^- 
spectively, (lA. 14p is changed by 



1 \ _V-. + V+ 

J ~ 2 ■ 



Hence, we have now change i/2 by 1/2 in the first equation of (lA. ISP and —i/2 by 1/2 in the 
second one. Respectively, (lA. 16^ for the second column reads 



B- 



B- 



Ak+ ' " Ak_ 

Then the second column Rjj{X,x,y) of the resolvent reads: 

Rii{X, X, y) = Tii{X, X, y) + P//(A, x, y), 

where 



Tn{X,x,y) 



ik+\x-y\ _ Jk+{\x\ + \y\)\ _ ik.\x-y\ _ ^fc_{|a;|+|^/|)^ 



Ak. 



- —;—[e 

Ak. 



and 



Pii{X,x,y) = — 



Note, that if ?/ < we get the same formulas. 



(A. 26) 



(A. 27) 



(A. 28) 
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A. 2 The poles of the resolvent 

The poles of the resolvent correspond to the roots of the determinant (lA. 20p . 

D{\) := + 2ia{k+ + k.) - Ak+k. - = 0. (A. 29) 

with k± as in (lA. 7p - (lA. Sp . Thus -D(A) is an analytic function on C\ C-_ U C+. Since there are 
two possible values for the square roots in k± there is a corresponding four-sheeted function 
D(A) analytic on a four sheeted cover of Cwhich is branched over C_ and C+. We call the sheet 
defined by flA. SP the physical sheet. 

We will reduce the equation flA. 29P to the solution of two successive quadratic equations. 
These can be solved explicitly but the process involves squaring and thus actually produces 
zeros of the function D{X) rather than of -D(A). Therefore we will then have to check whether 
or not the roots do actually lie on the physical sheet. 
Step i) 

Denote by cr = /;;+ + A;_ . Then 

(T^ = 2k+k_ - 2uj (A. 30) 

by ()A. 7p . hence flA. 29^ gives the first quadratic equation: 

o? + 2iao - 2{a^ + 2cu) - /J^ = 0. 

Rewrite it as 

2 o2 

a2 - laa = " ~^ - 2u =: 6 (A. 31) 

Finally, 

- = Y^f-T' (A. 32) 

where the root is choosen arbitrarily. 

Further let us express the roots in uj. Since a = 2y/uj, a = a + b/2, (3 = h/2 then substituting 
5 from ( lA. 3ip . we obtain 

8 = — -2uj= ^ '-^ = + — = 2a - hf < 0. 

442 4 82 4 16 4 16^^ 

Now flAT^ reads 

Q(Y 1 1 r -1 

a = -±-{2a-b) = - [(2a + b) ± (2a - b)\ = ^7,-, j = 1, 2, (A. 33) 

where 7^ e IR, and 

7i = a = a(C2), 72 = &/2 = a'(C^)Cl (A. 34) 



Step a) 

It remains to calculate the correponding spectral parameter A. First, the quadratic equation 
dAHni) implies by (lAT^ that 



On the other hand, 



A{k+k^Y = {2LU + = i2u; - 7J)', J = 1, 2. 
k+k- = V—oo + iXV~io — iX, 



(A. 35) 
(A. 36) 
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hence (lA. 35|) gives the second quadratic equation 



Therefore, 



A2 

Finally, we obtain four roots 



4 



A,=z|^4^-7|, (A. 37) 



where j G {1,2} and the square root can takes two opposite values. 

Corollary A.l The four-sheeted function D{X) has the following roots (zeros): 
i) j = 1 gives Ai = since 4a; = = 71. 

a) // I72I < 2^/cJ, then both j = 2 roots ±i|A2| are pure imaginary. 

Hi) If I72I > 2y/uJ, then both j = 2 roots ±|A2| are real: one positive and one negative. 



Remark A. 2 Note that a priori we can meet the wrong sign oflmk± squaring liA. 30) which 



is why the above calculation yields roots of D{X) rather than the physical branch D{X). Since the 
formulas ^A. 23\} - ^A. 281) involve only D{X) it is important to know which of these are actually 
roots of D{X) and also to know the multiplicities. This is done in the next two sections. 



A. 3 Discrete spectrum A = 

In order to check that the roots of -D(A) given in Corollary I A. II are actually roots of D{X) it 
suffices to check explicitly that D{X) vanishes (with the assumption that we are on the physical 
branch defined hj lmk± > for A G C \ C-t. 

For j = 1 we have 7 = 71 = a = 2^/^J and then Ai = 0. For j = 2 we have 7 = 72 = a'C^. 
If I72I = 2y/uj ( equivalently \a'\ = a/C"^) or 72 = ( equivalently a' = 0), we have A2 = 0. 

Let us check that A = is a root of -D(A): 

D{0) =a^ - + 2ia2i^ + 4cj = (a + 6/2)^ - h^j^ - 2{a + 6/2)a + = 
since k± = i^fuj. Now let us compute -D'(A): 

B\X) = ia{ \ + ~V ) - ( , . v/=Z7^+ ^^2L= ■ V-c^ + zA). 

\J —ijj ^ %X y—uo — iX y—uo + iX y—uo — iX 

Hence -D'(O) = and A = is the root of -D(A) of multiplicity at least 2. Further calculation 
shows that the Taylor series for D near zero takes the form: 

D[X) = {l-^)X' + 0{X% (A. 38) 

Therefore A = is the root of -D(A) of multiplicity 4 if and only if 6 = 4^^, i.e. a' = a/C"^, 
and we have proved the following lemma: 

Lemma A. 3 If a' = a/C"^ then X = is a root of the determinant D{X) with multiplicity 4, 
otherwise X = is a root of the determinant D{X) with multiplicity 2. 
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A. 4 Nonzero discrete spectrum 

Now let us check whether the roots A = A2 7^ corresponding 7 = 72 ^ {0, ±2a/cJ} he on the 
physical branch. We analyze two different cases: < I72I < 2^/0; and I72I > 
I. The case < I72I < 2^ (equivalently < \a'\ < a/C"^). 

Since 407 — 72 > 0> ^^e corresponding roots A2 are pure imaginary by (1A. 370 . Moreover, 
IA2I < cu. Indeed, (lA. 371) implies 

oj"^ - |A2|^ = + 7^/4 _ = (cj - 7^/2)2 > 0. 
Hence —uj =f iA2 < and k± are pure imaginary with nonnegative imaginary part, that is 

A;+A;_ < andlm(k+ + k_) > 0. (A. 39) 

The equations flA. 35^ and flA. 301) imply 

1 fl^ 
\k+k_\ = -\a^ -2{a'fC% {k+ + k_f = -2uj + 2k+k_ = -— + 2k+k_. (A. 40) 

In order to obtain and fc+ + k_ from the last two equations we have to divide the set 

< I a' I < a/C^ onto three subsets: 



a a 



{-a/C\ a/C^) \ {0} = {-a/C\ -a/^C^] U ((-a/v^C^, a/v^C^) \ {0}) U 

1) First consider the case a' G {-^ — , -^). Then ( lA. 391) and ( lA. 401) imply 

v 2C 

fc+A;_ = ^(a^ -2(a')'C^). 

{k^ + k.f = -^^ + ''^-{a'fC' = -{a'fC\ 
k+ + k_ = ia'C'^, 

and using (lA. 29|) . we obtain 

D(A2) = {a + a'C^f - {a'C^f + 2i{a + a'C2)(A;+ + A;_) - Ak+k. 

= + 2aa'C2 - 2(a + a'C2)a'C2 - + 2(a')'C74 = 0. 

Note that each 72 defines two values A2 up to factor ±1. If we replace A2 by — A2, fc+ and 
change places and our calculation remains valid. Therefore, both values of A2 are roots of -D(A). 

2) Further consider a' G (— t;;7t, ^ — 1. In this case 

k+k_ = ^{a^ - 2{a'fC^), k+ + k^ = -ia'C^. 

Then we have 

D{X2) =a^ + 2aa'C^ + 2(a + a'C^)a'C^ - + 2{a'fC* = 4a'C\a + a'C^) ^ 
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since a' 7^ and a' 7^ —a/C^. Therefore in this case both vahies of A2 are not the roots of -D(A). 
3) Finally consider < \a'\ < " ^ . Then OA. 391) - ( lATlOl imply that 

k+k_ = ~{a^ -2{a'fC^) < 0, 
(A;+ + A;_)2 = -a2 + (a')2c4<0, 



k+ + k-= i^a^ - ia'yC\ 

Then we have 

D{X2) = a{a + 2a'C^) - 2{a + a'C^)^ - {a' fC^ + a" - 2{a'fC\ (A. 41) 

To solve the equation -D(A2) = with respect to a', divide the RHS of ( lA. 411) by 7^ and 
denote p = a/C^ > 0. Then we get the equation 



p'^+pa - (a)^ = {p + a')^p'^ - {a'f, < \a\ < p/V2. (A. 42) 

Squaring both side of (1A. 421) . we get 

2(aY-/(a')' = 

The equation has no solutions for < \a'\ < p/V^ and hence D{\2) does not vanish. 
Corollary A.4 ^(Aa) = z/a' G [-^, ^). 
z.;Z)(A2)^0z/a'G(-^,-^)\{0}. 

II. The case I72I > 2y/uj (equivalently \a'\ > a/C"^). 

Since 4u; — 7I < 0, the corresponding roots ( lA. 371) are real: A2 < < A^, A^ = — A^. It is 
easy to prove that k± take the form: 

k± = ±/i + iz/, z/ > 0. 

Therefore 

= -/i^ - z/2 < 0, k+ + k_ = 2w (A. 43) 



1) First consider the case a' > a/C^.Then by flA. 40p and flA. 43^ 



k+k_ = ^{a^ -2{a'yC^), {k+ + k^f = -{a!fC\ k+ + k.=ia!C\ 

Therefore 

D{\2) = a{a + 2a C^) + 2i{a + aC^){k+ + k_) - Ak+k_ = 

a{a + 2a'C^) - 2(a + a'C^)a'C^ - + 2{a' fC^ = 

and then A2 are real roots of -D(A). Hence, the case a' > a/C^ is linearly unstable. 
2) Further consider the case a' < —a/C^. Then 

k+k_ = ^(a^ - 2{a'fC^) < 0, k+ + k_ = -ia'C\ 

D{X2) =a^ + 2aa'C^ + 2{a + a'C^)a'C^ - a' + 2(a')'C' = Aa'C^a + a'C^) ^ 
Therefore, in this case A2 are not roots of -D(A). 



31 



Corollary A. 5 i) In the unstable case a' > a/C"^: both A2 are roots of D{X). 
ii) If a' < —a/C"^ then neither of the A2 are roots of D{X). 

Summarizing, we have proved the following result 

Theorem A. 6 i) If a' E (— 00, a/(v^C^)) the only root of D{X) is \ = with multiplicity 2. 
ii) If a' G [a / \/2C'^ , a / C"^) , there are four roots of D{\): zero (multiplicity two) and ±z|A2| 
(pure imaginary) with A2 as in ( A. 31^ . 



Hi) If a' = +a/C'^, the only root of D{X) is X = multiplicity 4. 
iv) If a' G {a/C"^, +00), there are four roots of D{X): zero (multiplicity two) and ±|A2| with A2 
as in ( A. 37) . In particular there exists a positive root (linear instability) . 



Remark A. 7 Imagine reducing a' starting from a value greater than ajC'^ . Initially there are 
two real roots, ±|A2|, which approach zero as a' ^ a/C"^ from above, giving rise to an increase 
of the multiplicity of the A = root to four when a' = a/C"^ . As a' is reduced further below a/C"^ 
these two roots reappear as a pair of conjugate pure imaginary roots which move from zero to 
±iuj as a' goes from a/C^ to a/y^C"^. When a' = a/v^C^ these two roots touch the branch 
point (end of the continuous spectrum) and move onto an "unphysical" branch (on which the 
conditions M. ^] do not hold). As a' is reduced further these roots do not return to the physical 



branch and thus even when their magnitude becomes zero they do not coalesce with the physical 
A = root to increase its multiplicity and most importantly the spectrum is pure continuous 
apart from zero for a' < a/C"^ . 
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